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We study the spectrum of metric fluctuation in K-deformed inflationary universe. We write the 
theory of scalar metric fluctuations in the k— deformed Robertson- Walker space, which is represented 
as a non-local theory in the conventional Robertson- Walker space. One important consequence of 
the deformation is that the mode generation time is naturally determined by the structure of the 
k — deformation . 

We expand the non-local action in H 2 /k 2 , with H being the Hubble parameter and k the defor- 
mation parameter, and then compute the power spectra of scalar metric fluctuations both for the 
^vq , cases of exponential and power law inflations up to the first order in H 2 / k 2 . We show that the 

power spectra of the metric fluctuation have non-trivial corrections on the time dependence and on 
' the momentum dependence compared to the commutative space results. Especially for the power 

law inflation case, the power spectrum for UV modes is weakly blue shifted early in the inflation 
and its strength decreases in time. The power spectrum of far-IR modes has cutoff proportional to 
f**"» ' k 3 which may explain the low CMB quadrupole moment. 
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I. INTRODUCTION 



The history of the studies on the Cosmic Microwave Background (CMB) anisotropies and on the cosmological 
fluctuations is closely linked to that of the study of the standard cosmological model 0, Q. We now have high 
resolution maps of the anisotropies in the temperature of the cosmic microwave background y| , and its accuracy of 
the resolution is improving further. In relation to these observational data, overviews of the theory of cosmological 
perturbation applied to inflationary cosmology have been presented in Refs. The cosmological observations 

reveal that the Universe has non-random fluctuations on all scales smaller than the present Hubble radius. In the 
most currently studied models of the very early universe it is assumed that the perturbations originate from quantum 
vacuum fluctuations, which was first proposed in a paper by Sakharov Q. With this, the inflationary cosmology 
bears in it the 'trans-Planckian problem" : Since inflation has to last for long enough time to solve several problems 
of big-bang model and to provide a causal generation mechanism for CMB fluctuations, the corresponding physical 
wavelength of these fluctuations has to be smaller than the Planck length at the beginning of the inflation Q . Both of 
the theories of gravity and of matter break down at the trans-Planckian scale. Many methods have been proposed to 
cure the problem. The modification of the dispersion relation, which was used to study the thermal spectrum of black 
hole radiation Q , was applied to cosmology • Modifications of the evolution of cosmological fluctuations due to the 
string- motivated space-time uncertainty relations, <5x p hy S <^ > l 2 s: have been introduced by Brandenberger and Ho h0|. 
It was shown that the uncertainty relation plays a significant role in the spectrum of the metric fluctuation 
Greene et al. |l2j | proposed the initial states which give an oscillatory contribution to the primordial power spectrum 
of inflationary density perturbations. There have also been some attempts to explain the low CMB quadrupole 
moment contribution by using the pre-Big Bang scenario in string theory 01 . The ambiguity of the action in 
the presence of a minimal length cutoff in inflation by the boundary terms are studied by Ashoorioon, Kcmpf, and 
Mann 15]. However, it is not easy to construct a consistent field theoretic model which satisfies both the stringy 
space-time uncertainty relation and the spatial homogeneity and isotropy of the Robertson- Walker space. A direct 



'Electronic address: hckimOphya. yonsei. ac. kr 
^Electronic address: jhyeedphya. yonsei. ac.kr 
t Electronic address: rim@chonbuk.ac.kr 



2 



non-commutative deformation of the commutation relation, 

[x tt ,a?]=i8 1 "', (1) 

introduces a preferred direction in space, which breaks the isotropy of 3-dimensional space. Therefore, it would 
be interesting to construct a space-time non-commutative theory which keeps the spatial homogeneity and isotropy 
demanded by the Robertson- Walker space-time. 

Much attention has been given on the possibility of explaining the observational data as a quantum gravitational 
effect. As a theoretical framework to study these quantum gravity effects phenomenologically "Doubly Special Rel- 
ativity" (DSR, also called Deformed Special Relativity) was proposed by Amelino-Camelia where there exist 
two relativistically invariant scales, the speed of light and the Planck scale, and extensive studies have been fol- 
lowed [T7I ll8L fli^. I20L I21I I22] ] . Recently, it was argued that the coordinate space of the DSR theory defined in curved 
momentum space is described by the k— Minkowski space. Therefore, a good candidate to study the quantum gravity 
effect to cosmology is to extend the k— Minkowski theory to k— Robertson- Walker Space (k—RWS) and to study 
the effect of the deformation in the cosmological evolution. In n— Minkowski space the space-time coordinates are 
non-commuting generators of a quantized Minkowski space-time. The K-deformed Minkowski space-time introduces 
a dimensionful quantum deformation parameter, n, which can be chosen to have the dimension of mass |2.lL l24l l2^j . 
A natural choice of this deformation parameter is the Planck mass k = Mp. It is therefore important to construct 
consistent quantum field theoretic framework in k— Minkowski space, and explore the physical effects |2*oj in cosmo- 
logical evolution. In this respect, Kowalski-Glikman 26] have studied the effects on the density fluctuations of the 
quantum k— Poincare algebra. 

In Sec. II, we construct the theory of cosmological fluctuations in K-Dcformed Inflationary Universe (k-DIU). 
Starting from the scalar-gravity theory in a flat Robertson- Walker space we briefly summarize the linearized theory 
of scalar metric fluctuations. After developing the k-RWS, we write the theory of scalar metric fluctuations in the 
deformed space. We show that the scalar theory in k-RWS space is described by a nonlocal field theory in the 
conventional Robertson- Walker space. The nonlocal action is series expanded in H 2 /k 2 in Sec. Ill, where H is the 
Hubble parameter, and is quantized. In Sec. IV and V, we calculate the power spectra of scalar metric fluctuations 
for the cases of exponential and power law inflations, respectively. We show that the k— deformation alters both the 
time dependence and the frequency dependence of the power spectrum nontrivially. In Sec. V, we summarize our 
results. 



II. DENSITY FLUCTUATION IN k-DEFORMED ROBERTSON- WALKER SPACE-TIME 



If the Universe is quantum mechanically created with vacuum energy dominance, it will inflate from the beginning. 
Even though there may be other choices for the pre-inflationary universes, we assume that the inflation starts from 
the beginning of the universe. 

The calculations in this paper are carried out in 4-dimensional spatially Hat Robertson- Walker metric, 



ds 2 = -dt 2 + a 2 (t)(dr 2 + r 2 dn 2 ), 



(2) 



where a(t) denotes the scale factor of expanding universe. For later use, we introduce the conformal time rj defined 

by 

(3) 



»(*)■ 



The Einstein-Hilbert action for gravity coupled to scalar matter field is 



d 4 x. 



R 



1 



16nG 



dptp&>tp - V(ip) 



(4) 



where R is the Ricci curvature scalar. Starting from the action Q, it was shown that the scalar and tensor parts of 
the linear metric fluctuation are described by the action (See Refs. [HUlflj). 



5' 



di] 



(d v v k ) 2 + 



d 2 z 



(5) 



where / = 



d 3 k 



and 



H • 



(6) 
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with (fo being the scalar zero mode and H = a/a the Hubble parameter. It is noted that in cases of power law 
inflation and of slow roll inflation, H is proportional to tpo, hence z oc a. The action (jSJ) can be cast into the covariant 
form 

S = -\ j (Px^g» u dJld v K, (7) 

where 1Z is the gauge invariant metric fluctuation and plays the role of a massless scalar field in the inflating universe. 
The Fourier mode of 1Z, 7?-k, is related to the field by = zlZ^. Therefore, its power spectrum is 

^ - ^«°>- © 

We now want to write the action (JJJ in k— RWS. For this purpose, we develop scalar field theory in k-RWS in the 
following subsections. The 4-dimensional field theory in n— Minkowski space has been constructed in Ref. [rH 12(1 123| 
and references therein. We briefly summarize the results obtained in these references in the next subsection. 

A. «:-deformed Minikowski space-time 

The K-deformcd Hopf algebra H x describing the k— Minkowski space is generated by the coordinates x^ determined 
by the following relations: 

i 

[Xo,Xi] = -Xi, [Xi,Xj]=0, A(Xa) = Xu ® 1 + 1 ® Xu- (9) 
K 

The dual Hopf algebra Hk of functions on k— deformed four-momenta is described by the Hopf subalgebra of the 
k~ deformed Poincare algebra as follows: 

[k^kv] = 0, A(fci) = h ® e- ko ' K + 1 ® k h A(fco) = fc ® I + f ® k . (10) 

The first Casimir operator of the algebra (jHJl and (|f OJI is 



M 2 = ^Ksinh^J -k 2 e fc °/ K . (II) 

It follows from this that for positive k the on-shell three-momentum is bounded from above by 

k 2 < k 2 , (12) 

and the maximal value of momentum results in an infinite energy |27| . 

By using k— deformed Fourier transform, the fields on k— Minkowski space with non-commutative space-time coor- 
dinates x — (xq, Xi) is written as (kx = kiXi — k x a ): 

$(i) = J k ®(k) ■ e lk± :, (13) 

f f d A k 

where = -, — rj and $(fc) is a classical function on commuting four-momentum space k = (ko, fcj) and the normal 

Jk J ( 27r ) 
ordering is defined by 

• e ikS: '= e ~ ifc oio e ikx ^4) 

Multiplication of two normal ordered n— deformed exponentials follows from Eqs. © and (|14fl : 

■ e ikS: " e iqi ■=■ e i ( ke_, ° ' K +q)x-i(feo+?o)io . 0-5) 

which follows from the four momentum addition rule described by the coproduct I|1U|I . From this we get the conjugate 
field, 
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$ f (i) = / &(k) : e ik£ :, $ f (k, ko) = e 3k °/ K $*(-e ko/K k, -ko). 

Jk 



(16) 
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The multiplication of fields can now be expressed as 

<P 2 (x)d 4 x = / i>(fc)$(-ke fc °/", -fc ). (17) 



The differential calculus and its covariance properties under the action of k— deformed Poincare group have been 
constructed in Ref . [2^, |2f| . The left or right partial derivatives &a to define K-deformed vector field are given by 

d A <$>(x) =: X A (j£U $ (x) :, (18) 

where \A '• e%kx XAi^^e 1 ^ : and 

1 1 2 

Xi{K) = e Wre fci, Xo = «sinh^ + — e k °' K . (19) 

k 2k 

The adjoint derivative d\ can be defined to satisfy 

' d 4 i$!(x)5 $ 2 (£) = / d 4 i[^$ 1 (i)]$ 2 (£), (20) 



which leads to 

dt$(x) = xlidJiMx), (21) 

where 

xi(fc)-XM(-e fco/K k,-fc ). (22) 

Based on these, the A</> 4 field theory was constructed in Ref. 23]. In the next subsection we generalize the formulation 
to the case of the Robertson- Walker space-time. 



B. The scalar field theory in K-deformed Robertson Walker space-time 

To generalize the field theories in the K-Minkowski space to the curved space case, we must be careful in selecting 
the coordinates which satisfy the commutation relation Since any non-decreasing reparametrization of t is an 
equally good time coordinate in commutative space-time, it is important to choose the time coordinate for which the 
commutation relation, 

[xo,Xi] = -Xi, (23) 

K 

is imposed. We note that a natural time coordinate consistent with the commutation relation is the cosmological time 
xq = t. This choice ensures the same form of commutation relation satisfied by the locally flat coordinates (i, Xi — 
a{t)xi): = a(t)[t,Xi] — iXi/n. This simplicity cannot be attained for other choices of time coordinate. For 

example, consider a time coordinate xo defined by the 00-part of the metric goo = —s 2 (xq). In the locally flat 
coordinates, the commutator [s(xq)xo, a(xo)xi] — a(io)s(io)[^o, Xi) + a(xo)[s(xo),Xi]xo is not simply reduced to a 
well defined form of Eq. 10 . In this sense the natural choice for the time coordinate is the cosmological time t where 
s(i) = 1. With the cosmological time t, the equations (P|)^ i|18|) can be used without modification. 

The generalization of k— deformed vector fields in k-RWS can be written as Eq. I|18fl with the operator x^ defined 
by, 

Xi = e d °/ iK ^, xo = Ksinh — - —d t e da ' 2lK g l3 d 1 e da ' 2lK , (24) 
i in Ik 

where the ordering of the time-dependent metric, g lJ , and derivatives are determined by demanding the adjoint 
derivatives to satisfy Eq. I|2U|I . which gives 



x\(k) - -k h (25) 
ik v ' 2k 



dl<f>(x) = -Ksinh— - ^-e do ^ 2i ^di\e do/{2iK) d 3 <f>(x)\ g l \t). (26) 
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In addition to Eq. (|20|l we demand the condition 



(27) 



to determine <9g. 

Given the covariant derivatives and its adjoint derivatives, we can write the action of a massless scalar field in 
k-RWS as 



(dl&(x))a 3 {t)d $(x) ~ (SI&{x))a{t)di<Z>(x) 



(28) 



where we choose the symmetric form in the action so that the metric dependent factor is placed in the middle of the 
operator products. It turns out that this choice gives the simplest form of the interaction between different modes. In 
addition we do not consider the change of measure psl ] due to the complication of the non-commutative multiplication 
since what we are interested in in this paper is to understand the main feature of the k— deformation on the metric 
fluctuation. 

Using Eqs. (O, JT^), lO, and we obtain for the action 



S = ^ / *<{ 5 (k,t)$_ k (t-^)|. k (t+^)- a 3 (t)[p 3 (t)+ 3 3W^W]*-kW*kW 



(29) 



where k denotes the relative ratio between the physical momentum (|k|/a) and the non-commutative scale (k(t) = 



i(t)n 



) and the coefficients g(k, t), gi(t) and p n {t) are given by 



ff(M) 
31® 
92 (*) 

93(t) 
Pn(t) 



a 3 (t) 

1 
2 



l + (.gi(t)-2)fc 2 + 



a 2 {t) 



+ 



92® 



(30) 



a 2 (t + i/(2 K )) a 2 (t-i/(2n)) 



a 2 (t + i/(2K))a 2 (t-i/(2n)) : 



a 3 (t — ijn) 



a 3 (t + i/k) 



a(t)a 2 {t-i/(2n)) a(t)a 2 (t + i/{2n)) 
a n (t-i/n) a n (t + i/n) 



j»(t) 



a n (t) 



Note that the left hand sides of Eq. I|30() are defined to satisfy lim <?; = 1 = lim p n . The action (|29() is highly 

Kt — >00 Kt — >O0 

non-local in that the fields are non-locally multiplied in the action, in addition to the nonlocal coupling between the 
background metric and the field modes. Each mode of the field, $±k, is diagonalized so that it is not coupled to other 
modes of different k. 



III. FIRST ORDER APPROXIMATION AND THE HAMILTONIAN FORMULATION 

It is not possible to solve the nonlocal equation of motion derived from the action (|29|l exactly. The canonical 
formalism for Lagrangians with non-locality of finite extent has been proposed by Woodard |29j. However, we do 
not follow the formalism since our purpose is to obtain information on how the non-commutativity (|23|l affects the 
evolution of the metric fluctuation in inflationary Universe in a simple calculable form. Instead, we use a perturbative 
expansion in the parameter H 2 = H 2 /k 2 and construct the Hamiltonian for the action up to the first order in H 2 . 

To have an approximation of the action l|29|l , we expand the integrand in H 2 as 



g(t)<Z>(t+-)<Z>(t-~) 

K K 



(-3) 



-,(-2) 



6(^ 4) 

_ 9 
4k 4 



$ 2 (t) + 0(H 6 



(31) 
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where g( n '(t) denotes the n th indefinite integrals of g(t), and 



g ( s n) (t) = - (t + %/k) + 9 ( - n) (t - i/K)] , g [ A n) (*) = ^ [9 { ' n) (* + */«) - .9 ( ~" } (* - */«) 



From Eqs (|29[1 and l|31|l. we get the action up to the order 0(H 2 



S 



dt a d (t) i /x(k, t)$ k (t)$_ k (t) 



.94 



55 fc 5 



where the coefficients are given by 



Pi + 7 1 .93 - 



a 3 (i + i/K) 



a 3 (i - 



(32) 



(33) 



(34) 



a(t)a 2 (t + 3«7(2k)) a(i)a 2 (t - 3i/(2«)) 



.95 (t) 
7(M) 



a{t)a 3 (t + i/K) 



a(t)a 3 (t — 



2a 3 (t) 



a 2 (i + i/2 K )a 2 (t + 3i/(2«)) a 2 (t - i/2fi;)a 2 (t - 3z/(2k)) 
_ 6/t 4 (9 (-4)_^-4) ) + K : 



The asymptotic values of these coefficient functions are 

54(00) = 55(00) = fj,(0, 00) = 7(0, 00) = 1, 



(35) 



which make it easier to guess the asymptotic behaviors of the coefficient functions for large t. 

For notational simplicity, we use the change of variables $k,+ = gC^k + $-k), ^k- = §(3?k — $-k), to write the 
action in a diagonal form: 



S = 



1 



dt a 3 (t) M (k,t)$2( t )_ 



.94 



55 k 2 



2 7(M$3 
Q 3k 2 q 



(36) 



where a = (k, ±). Note that the coefficients fx, 7, and 5, are exactly calculable once a(t) is given. 

Introducing the conformal time r] and the rescaling of the field, <p a — a(t(i]))Q a , we reduce the action (|36fl into the 
form, up to the order H 2 , 



S 



- / dr 1 ifi(d v <j> a y-u(k,t)<Pi + -^ 



7CM) 



3K 2 a 2 (i) 



where 



AjH 2 



5H 



3n 



3k 2 



AH 2 4#7 



w(k, t) = fc 2 f 54 

A 2 = a 2 (77) 



5 5 P\ 5 r) (A 2 ff) , 



4 

/i(k,t) 



a 

d v ('~fdr l a) 



3n 2 a 



(37) 

(38) 
(39) 



3« 2 a 3 



Note that the action 13 7|) contains a higher derivative term, which may lead to nonunitary evolution of the system. 
Since this higher derivative term is a term of order H 2 and our purpose in this paper is to obtain the effect of the 
deformation on the cosmological evolution up to the I st order in H 2 , we require that the higher order derivative term 
is written as a function of the field, its first time derivative, and time: 



a 2 ^ = M>(0,^0,?7). 



Explicitly, we use the linearized approximation, 



a{rj) 



7 1/2 (M) 



[c(v)<f> + d(r))dn<f>], 



(40) 



(41) 
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where the coefficients c and d are to be determined by consistency. 

This requirement is equivalent to the perturbative calculation up to the 1 st order in H 2 . This can be shown as 
follows: The equation of motion for <f> — (f>o + -^s4>\ + ■ ■ ■ can be written as 

d v {fLd n (/)o) + lo^o = 0, (42) 

dr,(fid v <l>) +ucf> = ^d 2 (^<9 2 o ) , 

where the first equation is the th order equation and the second is the full equation written explicitly up to the 1 st 
order. The first equation of Eq. 1(42(1 defines d 2 (f>o as a linear function of d v 4>o and 4>q. Then, the second line can be 
understood as a defining equation of d 2 (f> as a linear function of d^fy and <f> up to 0(H 2 ), which is Eq. 141(1. 

With this reasoning and Eq. 1)41(1 . the action 1(37(1 is perturbatively equivalent up to H 2 to the following unitary 
action 

S=\j dq [m{d^ a f - f(k, rfifc] ) (43) 



where k = Ikl and we have 



m = ^ + ( 44 ) 

Substituting Eq. ((41(1 into 1(3 7|) and requiring the resultant action to be the same as the action ((43(1 with ((41(1 as its 
equation of motion, we find 

/7 l/2 W ,2 

c= . a— . (45) 

ma to 

Eqs. ((4411 and l(45|) imply that m and / satisfy the following differential equations: 



(46) 



7 (to) 2 
3k 2 to 

2—5-2 + / = w ( k > *) H 



3k Tn a 



Note that these conditions make the action to be a functional of d v 4> and 0. The time evolution for the theory l(43|) 
is unitary and quantum mechanically well defined. 
We introduce mode dependent conformal time rjk by 

drj k = mT l (k,ri)dr) = ^ . (47) 
a(t)m{k,t) 

Then the action 1(431) can be written in a simplified form: 

S=lj d^^-nUrjk)^}, (48) 

where ' denotes the derivative with respect to rjk and 

nK»fe) = m(M fc )/(Mk)- (49) 
The Hamiltonian for mode a is the same as that of the time-dependent harmonic oscillator with frequency squared 

ff« = ^f + ifi2fa*)& (50) 
The time evolution of each mode can be described by introducing invariant creation and annihilation operators > 

^ a = ~^/2^f 3 *cfi a ~ fa'M, Al, = -^((p a n a - ip a '$ a ), (51) 



8 



where ip a is the mode solution of the differential equation (|54() below and A a and A^ a satisfy the Liouville-von Neumann 
equation, 

m r , h A a + [A a ,H a }=0. (52) 
One may invert Eq. (|51|l to construct the field operator in terms of the creation and annihilation operators as 



h 1 ' 2 
h 1 ' 2 



tp a {-nk)A a + ^*(%)i^ 

<p' a (Vk)A a + ip*J(r)k)A\ 



(53) 



Note also that the Liouville-von Neumann equation is equivalent to the following differential equation for the coeffi- 
cients (fi a , 



The commutation relation [A a: A^] = 5 a p restricts the mode solution ip a to satisfy fa^a' ~ P'aVa = 
We present the first order approximation of m and / for later use. To first order in 1/k 2 , Eq. (|46|l gives 



(54) 



■^2 



TYL — fl - 



7 (m) 
3k 2 p 2 



(55) 



/ = MM 



1+ 1 



4uj 



3n 2 m 2 a 2 



1/2 



w 1 



fl- u ) 
V 3K 2 m 2 a 2 1 



Using the explicit form for w and fl, we have 



H 2 



m ~ fi ~ 1 - -^-(1 - 7ei) - fc 2 , 



(56) 



where e r , 



1 + 5**' 



2i7 2 a 



2„2 



ei # 2 



2 , £ 2 



1+ 2"""6~ (1_16ei_£ ' + ^ +e3 



with 



(ff)n+l ; 
inflation. From these we have 



-, are constant numbers for power law inflation and vanish for exponential 



-in s 1 ff 2 fc 2 
m (fe,?7) ~ 1 + a n — ~- + —o—n H , 

ilUrjk) ~ fc 2 (l-«; 1 i7 2 - U ; 2 P)-(2 + e 1 )i/V(r / )(l-u )3 5 2 ) + 
For the power law inflation, the values of a n and Wi are given by 



(57) 



<h, = —zr—-, w\ = i(13 - ll/n), W2 = 19/12, w 3 = 



1 



If 45 7 3 



6n ' * 6 V " '"" * "'" ' ° 3(l-l/(2n))V 4n An 2 n 3 
and for the exponential inflation their values are given by the limits n — > oo of Eq. (|58|l . 



(58) 



IV. METRIC FLUCTUATIONS IN k-DIU: THE EXPONENTIAL INFLATION 

The simplest inflationary model is the exponential inflation, in which the scale factor a(t) increases as, 



a(t) = ape 



HI 



OO < t < oo. 



Here ao is the scale factor at t = and H is the Hubble constant. Using the conformal time 77, we get 

1 



Ht = — ln(— agHri), a(rj) 



-Hr) 1 



(59) 



(60) 
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where the conformal time r\ varies from — oo to as t varies from — oo to oo. From Eqs. (|34(l and l|38(l . we have 

sin if 2 



sin3if 4 - 9 , — . 



3H 

7 (k,i) = £(3fl)-£(5) 



cos if 



A' 2 



p 

cos if — — 



(61) 



where £(x) = 3 



1 — cos x sin x 4 cos x + 1 1 
-3 h 4 

x A x 3 2x 2 8 



13 

80" 



Since the action H48|l is obtained in the ff 2 expansion, the normalized Hubble constant, H = H/k is assumed to be 
smaller than one. Moreover, the Eq. I|67|l below restricts k = k/{a{t)n) to be smaller than one. Then we get m(k,r]) 
from the differential equation (|46|l and (|61|) by series expansion in H 2 and k, 



H 2 jx 2 
3 if 2 /2 2 



(1 - if 2 fcV ) + 0(ff 4 ). 



By integrating m 1 over r\ using Eq. I|47(l . we get the mode-dependent conformal time rjk, 



In 



Hkr]\ 



2Hk(l-H 2 /6) \l-Hkr)J l-H 2 /6 



1 



if 2 



(62) 



(63) 



which is normalized so that rjk — at r\ = 0. A crucial point is that there is a global rescaling of the conformal time 
due to the non-commutative effect. 

From Eq. (|57|) . the effective frequency squared il 2 is written as 



VliVk) = k 2 (l- Wl H 2 -w 2 k 2 )-2H 2 a 2 (r 1 ){l-w 3 H 2 



(64) 



1/4 



w 2 E 2 kSl 



4 



where W\ = 13/6, u>2 = 19/12, w 3 = 1/3, and 



k 2 = k 2 [1 - (wi + 2/3)H 2 ] , v 2 — 1/4 



2(1 -w 3 H 2 ) 
(l-iJ 2 /6) 2 



2. 



(65) 



The explicit value of v is v = 3/2 to this order, which is the same as the commutative space result. We note that 
both the frequency and the mode solution have corrections from non-commutativity for large |?7fe|. The last term of 
the second line in Eq. (|64J) becomes negligible for rjk ~ 0. 



A. Mode generation and the initial condition 

In the k-RWS, spatial momentum is also restricted similarly as in l|12|l . In terms of comoving momentum k = |k|, 
we have 

4t < «• (66) 
a(t) 

This gives the upper bound of k 

k < k max (t) = Ka(t). (67) 

The maximal value of the wave-number l|67ll is very similar to that used by Brandenberger and Ho ^(| except for the 
fact that the maximal value used in Ref. [T(| is determined by an effective scale factor modified by the Moyal star 
product in the action. 

Since m is positive definite, the mode-dependent conformal time rjk, Eq. I|47l) . is well defined and is an increasing 
function of t. The relation l|67|) implies that for a given fc, there exists a conformal time 77°, the time saturating the 
relation l|67|) : 

a(t(4)) ^ - (68) 

K 
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Then, the mode 4> k cannot exist before 77^. In other words, 77° is the generating time of the mode tp k . This provides a 
hint to one of the major issues in which state the fluctuations are generated. To satisfy the continuity of the number 
of quanta of the k mode when the mode becomes physical at 77°, it must be in the adiabatic vacuum state. This 
vacuum state can be chosen to be the WKB mode solution, 

lira ifkim) = , 1 , ; expi ( / Slkdr) k +ipk ] , (69) 
m^v° k V2fi fe (?7fe) \J n ° J 

where the constant phase ip k can be chosen conveniently 



B. inflationary evolution 



The inflationary evolution of the mode solution is determined by identifying f2| of Eq. I|64|) . The corresponding 
commutative space values can be obtained by setting H = in Eq. (j64(l . With the f2| we have two different time 
scales and jf k defined by Q' k (ri k ) = and Q k (ril) = 0, respectively. These time scales are given by 

_c_ ~ V4) 1/4 < - 1/4)V3 , (y8 _ 1/A)W2S 2 x 

Vk ~ ~ fa^iTV ' Vk ~ k \ 2 + ) ' (?0) 

where Q k (i] k ) increases while r\ k < r] k and decreases later as shown in Fig. 1. It is positive definite when r/ k < rfL 
and negative later. For r/k > rj k , the modulus of 4> k increases in time. By Eq. 1)68(1 . the mode 4> k is generated at 
77° ~ —[kH)~ l . Note that these three time scales satisfy 

^ < ^ < Vk, (71) 
if H <C 1. Note also that during 77° < 77^ < rfj., the condition for the WKB approximation, 

^^u^Mi^l, (72) 

^ l k 

is valid. In Fig.^ we present schematic plot of fl k (j] k ) for a given mode <p k . Therefore, during this period, the WKB 

«2 




FIG. 1: Schematic plot of Qj? for the exponential inflation and for the UV modes in the power law inflation. There is no mode 
4>k for rjk < Vk- The mode solution for r)1 < r\k < rj^. is given by the WKB solution and the mode solution for 77^ > rj% is given 
by the Bessel functions. The two solutions are matched at 77J:. 



mode solution, 



fVk \ 

ifik(Vk) = ■== expi I / tt k dr) k + ip k , (73) 
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can be used to describe the time evolution. We use this solution to determine the matching condition at r] k = 

^,) = -i=, v'M)=i\[^, (74) 



where f2 2 = ^(Vk) — k 2 [\ — 2(y 2 — 1/ 4) X / 2 H], k given in (68), and tpk is chosen to give this matching condition 174(1 . 

For r\ k > r\ kl the last term in Eq. (|64[) becomes much smaller than other terms. Thus we ignore this term and use 
the Bessel function as the solution for rj k > i] k , 



Matching the two solutions at r\ k = rj^, we get 



A k = - 



B, 



z v Vk 



(75) 



(76) 



Note that A k and B k are independent of k since kv^ is independent of k due to Eq. (|70|l and f2 c oc k. 



As kr] k 



0, the Bcsscl functions become 



Jv 



{-kVkT 
2 v v\ 



2 v {v- 1)! 
ir{-k-q k y ' 



(77) 



and the second term in Eq. I|75|l dominates in the later time (—kr/ k ~ 0). Therefore, the power spectrum of the scalar 
metric perturbation has the form 



Pk{t) 



k 3 \cp k (t)\ 2 [2>-l)!] 2 |f?| 2 tf 2 (l + H 2 /6f 



2vr 2 z 2 (t) 



27r 4 (z/a) 2 



(1 - 17i? 2 /6) 3/2 



krjf 



3-2v 



(78) 



where z/a, H, v, and B are constant numbers. Note that 3 — 2v ~ §(103 — 2a)H 2 = in the present case since 
u>3 = 1/3. Therefore, the spectrum of the metric fluctuation for exponential inflation in n— RWS is time independent 
and is scale invariant up to the first order in H 2 . The only effect of the non-commutativity to the power spectrum 
is the global rescaling of the power spectrum, which is of the order H 2 . It is an interesting fact that W3 and 2a 
are the same. Note that a originates from the scale factor of the mode-dependent conformal time r] k with respect 
to the conformal time 77, and W3 comes from the H 2 order correction term of the frequency squared. Since there is 
no physical reason for the coincidence, it is possible that the next order correction may give a result of W3 > 1/3. 
This is an interesting possibility since this positive power of kr\ makes the power spectrum decrease in time. If the 
present analysis is applied to the tensor mode fluctuation, the time dependence can be used to solve the gravitational 
hierarchy problem (H/Mp ~ 10 ) 31] . This is what happens in the power law inflation considered in the next 
section. The spectrum (I78f) is scale invariant in contrast to that of Ref. 26] with the same initial vacuum state. The 
difference may be attributed to the different choice of the initial conditions. At the present case, the initial time is 
dependent on the mode through Eq. (|68|l . which uniquely fixes the initial state l|69|l for the mode solutions. 



V. METRIC FLUCTUATIONS IN k-DIU: POWER LAW INFLATION 



In this section we calculate the metric fluctuation in the power law inflationary model, in which the scale factor 
increases as, 

a(t) = a (Kty\ < t < t f , (79) 

where n ^ 1, ao is the scale factor at the Planck time t = 1/k, and tf is the instance when the inflation ends. In this 

[2 

model, the variable z(t) in Eq. JBl is given by z{t) = \ —Mpait). For n 7^ 1, we have 

V n 



Kt = I — 

,V0 



a>(v) = 00 I — 



(80) 
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3n — 1 

where fi = — — and r/o is the conformal time corresponding to Kt = 1, given by 



2(n-l) : 



?7o = 



^ - 3/2 
a n 



(81) 



*7/ = To( K */) 1 ^ ( - 2a ' ^ is the time when the inflation ends. 
When /cf ^> 1 is large and k 2 /(a 2 n 2 ) <C 1 is small, we have 



L (fc,7y) ~ 1 + On 



a 2 (rf)n 2 



The mode-dependent conformal time r\ k can be approximated as 



Vk-V 



2{p-l)' 



=2/x-3 



where we use the notation 



V 



V 



2/iK 2 aQ 



7.2fJ.-l 



(82) 



(83) 



(84) 



rj 

and we normalize the time so that rjk = when r\ = 0. The function £l 2 in Eq. (|57|l for H 2 <C 1 is approximated to be 



n 2 k {r, k ) ^ k 2 {l-w 1 H 2 -w 2 k 2 )-^-^H 2 a 2 (r ] ){l-w 3 h 2 ) + 5, 



(85) 



where S represents the smaller terms proportional to the differences of the Hubble parameter from its time-averaged 
values, 

S = - Wl k 2 {H 2 - Hi) + 2w 3 H 2 a 2 (H 2 - h 2 ). (86) 
WiS are given in Eq. (57) and the time-averaged values of the Hubble parameters, Hk and h are defined by 



Hi. 



Jdy k H 2 ( v ) 
k 2 J drj k 



h 2 = I di] k a 2 (r))H 4 (r]) 
k 2 J drjka 2 (r])H 2 (r]) 



(87) 



with the rjk integrations performed over the validity range of the differential equation (|54|l for a given mode solution, 
which will be clarified in the next subsections. We do not put the index k to h since h depends on k very weakly. 



A. Mode generation and the initial condition 

As in the case of the n— deformed exponentially inflating universe, due to the condition (JB7J, the mode 4>k is 
generated at the conformal time rj(k), 

(l. \ ~ f»-l/2 
— , (88) 

a nj 

where we assume the Universe is inflating with n 1. 

A serious obstacle in finding physics of low comoving momentum modes is that the action (I48|l is not well defined 
for large H 2 since the action is approximated by expansion in H 2 . The condition H(rj m ) ~ 1 is attained at 

Vm ~ Vo . (89) 

Before ry m , our approximation for the action l|48|l is not valid. This condition restricts the validity range of the present 
approximation to the modes (f>k with comoving momentum 

k > k m = n^- 3 / 2 a Q K. (90) 

In this subsection, we restrict ourselves to the ultra-violet modes satisfying the condition l|9(J|) . To know the behavior 
of smaller frequency modes, one should use better approximation of the action (|29|l instead of (|48p. We present some 
reasonable arguments for the evolution of those low frequency modes in Sec. V.C. 
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For modes k > k m , all the arguments for the initial state in Sec. IV. A hold true. Therefore, the initial state is 
given by the WKB ground state, 



(fikivk) = , 1 , = expi ( / Clk(r)k)dr)k +ipk], k> k m , 



where r\ k is close to the mode generation time rj® given by Eq. I)68[l. 

*(*) 



\ 


\ \ 




\ 

\k 


\ \ 

\ ki N k 2 
\ 




1 FIR V sIR \UV 
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\ 






\ perturbatively unreliable 


pei'tsarbatively reliable 















4 



V kl Vn 



'lk 2 



(91) 



FIG. 2: Schematic plot for the mode generation time in power law inflation. Each curve describes the value of the normalized 
physical wave-number k(t) = — for a given conformal time -q. Since a(t(rf)) increases for expanding universe, the value of 
k always decreases. Due to Eq. 1881 each mode becomes physical when it crosses the horizontal line k — 1. We denote the 
unphysical part of each modes using the dashed curve. Since H = 1 at T] = rj m , the region r\ > r\ m is perturbatively reliable 
and the region 77 < rj m is perturbatively unreliable, which are divided by a shaded vertical line. The UV modes resides in the 
perturbatively reliable region. The IR and the Far Infra-Red ( < 1) modes pass through the perturbatively unreliable 
region. 



B. inflationary evolution 



The inflationary evolution of the mode solution is governed by fi|. When k > k m , we approximate fif. in Eq. I|85|) 

I the term 5 in Eq. I|85[) . as 

n k ( Vk ) ~ k - 2 -2-2% +••• • ( 92 ) 

tin 



using Eqs. (f5^f) and (JHSJl, and dropping the term S in Eq. JHSJ), as 

,,2 1 I A „„_ lA 



4 



where 



/ 2n — 1. o 

1 - (wi + )H 2 k 

[in 



H-l 



w 3 h 



(93) 



Eq. H92fl looks similar to Eq. I|64|) except for the final term. With the Vl\(rj k ) in Eq. (|92|l . we have two different time 
scales defined by Q' k {Vk) = and Q\(rf k ) = 0: 



V k = Vo 



v 2 - 1/4 



w 2 {n- l/2)( M -3/2) 2 fc4 



Vk 



^ 1/4) 



l + toa(Ai-3/2) 



1/(2^+1) 



v 2 - 1/4 
0-3/2 



(94) 



2/u 



p"2M 
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where fco = k/(nao). Since fco ^> 1 for modes under consideration, we have 



rf k « Vt « Vl 



(95) 



Since the condition for the WKB approximation, f^'/fi 3 . <C 1, holds during rjk < rfy we may use the WKB solu- 
tion l|91|) . Therefore, we get the mode solution at 77^ given by Ea. (|74|l with 



<? 2 



2n+ 1 v 2 - 1/4 



(96) 



where the second term in the parenthesis is smaller than 1 for fco ^> 1. For r/k > 77^, we ignore the last term in 
Eq. (|92[1 since fjk is very small there. The solution for the differential equation (|54|l is given by the Bessel function in 
Eq. I|75|) with parameters given in Eq. (|93|) . With the initial condition l|74|) and O c in Eq. Ij96|l . Ak and Bk are given 
by Eq. fty. 

Since we are considering modes with k ^> 1, we always have \kr]%,\ 3> 1. Using the asymptotic expansion of the 
Bessel functions, 



Mx) - v h cos[x - { " + Mx) - v h sin[x {v + 



we get 



At = ^./^£ e -«(fc^ + ^+l/2)V2) 

2 V k 



(97) 



(98) 



Note that the absolute values of and B(~ are very weakly dependent on k since f2 c 
From this, we have the power spectrum, 



A:. 



Puv,k{t) 



fc 3 \Mt)f 

2ir 2 z 2 {t) 



i[2"(v-l)\\ 2 (p-3/2) 1 - 2,/ / k 



M 2 



16tt 3 



QqK 



3-2v 



Since we are interested in the time evolution for 77^ < 77^ < 77/, we have the time-averaged Hubble parameters, 



(Vf) 



4-2(1 



{fj c } 4-2(i 



1 



2(2 - f i) 



-2(1-4 
If 



■Hf-Vk 

- m 2 ^ 



2(2 -n) nf-ini)- 1 " 2(2 — /J,) 



2(2 -n) {f\ c k Y^ 



(99) 

(100) 
(101) 



The mode-dependent conformal time 77^ at the end of the inflation is almost independent of the comoving momentum 
k. Note that from Eq. 



3-2i/ 



/1 — ch , c = 
3 - 2/i + 2ch 2 



- 1/4 



2// 



w 3 



(102) 



n- 1 



2c/i 2 



where c is a positive number for n > 1.53. 3 — changes sign from positive to negative at % = r)fi a t, where 

"(2- M )(2 M -3)" 



Vflat 



'/()■ 



(103) 



Since |3— 2u\ -C 1 always, the slightly blue-shifted spectrum changes into the slightly red-shifted spectrum at 77 = r)$i a t- 
The power spectrum decreases in time since /i — v is positive even though it is very small. For example, for n = 14, 

k 2 

we have a(?7/zat)/ a o ~ e an d -ffc ~ 0.075 — j. Since the currently discussed minimal duration of inflation is roughly 

Alp 

a/ ao ~ e 60 ?7/ must be later than r)fi at and the size of the spectrum decreases in time. We will observe the 
decreasing red spectrum of the scalar density fluctuation for high comoving momentum k > k m . 
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If n = 19, on the other hand, we have a(r)fi at )/ao ~ e 60 and we will observe fully scale invariant power spectrum 
for the high comoving momentum modes for minimal inflation. 

There are a couple of features in our result 199JI which can be used to resolve the gravitational hierarchy problem, 
the extreme weakness of the gravitational wave (the tensor mode fluctuation) relative to the large scale density 
fluctuations. The first is the fact that the power spectrum decreases as k increases for later time, which is also the 
case in the commutative space results. The second is the time-dependence of the power spectrum, which is a new 
feature of the present paper, that may describe the difference of the power spectra between the structure formation 
time and the observation time. These features may solve the gravitational hierarchy problem without imposing any 
fine tuning on H or on the non-commutativity scale k. 



C. Speculation for IR modes 



In the previous subsection, we deferred the discussion on the evolution of the modes of smaller frequencies than 
k m . Due to the definition of k m , l|90[) . these modes always satisfy in its physical region, 

<1. (104) 



a(t)n 

The Far Infra-Red (FIR) modes are created at rj = 77° < r/o. For 77^ < 77 < r) m , the condition 11041) can be interpreted 
as 

k 2 < a 2 H 2 , for 7/ < rj m , (105) 

since rj m is the time for H/n ~ 1. 

To figure out the dynamics of these modes, we must have reasonable approximation on the action around H 2 /k 2 ~ 1. 
We assume that the action l)48|l is still valid even though the explicit functional forms of m and / are not known. We 
keep Eq. (|46[) since the equation comes from the requirements of quantization and unitarity of the time evolution. To 
show that the mode solution at time rj m is almost independent of fc, we use several steps of reasoning. First, we show 
that m(k,r]k) is very small until 77^. < r\ m . Next, we argue that fl 2 is also very small. Finally, we argue that f2| is 
almost independent of k for this region of time. With these, we have the result that the matching condition of ipk(rik) 
at r]k ~ 77 m is almost independent of k. 

m(k,r]) is defined by the first order differential equation (|46|l with to = 1 at nt — > 00. The value of to increases 
in time asymptotically approaching to, to ~ ft at 77 = 0. Therefore, to get the behavior of m, we need to analyze its 
behavior around to ~ 0. In this case, the differential equation (|46[1 becomes 



± (106) 
to y 7 

Thus, to have a consistent, real-numbered value of to, we must have /i/7 < in this region and 

-1 



m{k,rik) 




A - / W drj k 



(107) 



where the integration constant A is a large number which ensures the value of to to be small. The value of to increases 
to O(l) around /i/7 ~ since for positive \x we have to ~ jl as a th order approximation. Let 77^ be the value of r]k 
satisfying P>(ti^.)/'y(t]^) = 0. Then, we have to ~ 1/A -C 1 for 77^. -C ??£ and to ~ p, for rjk > 77^. Rough estimation on 
the value of 77^ is possible by using the asymptotic form of m(k, r]k) ~ j2. Since 

l-On^-^<l, (108) 



K 2 a?{t c k ) 



where f? is the time corresponding to the mode-dependent conformal time 77^, and the mode 4>k is generated before 
77^. Since H(r/ m ) ~ k, we must have 77^ ~ r/ m + • ■ ■ . Therefore, the zeroth order of 77^ is independent of k. 

Next, we consider the function / of Eq. (46) for 77fc <C 77J:. From Eqs. I)46[l. H106[) . and the change of variable 

/ = —o /, we get 



(109) 
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where we have ignored the term proportional to m since m is very small. For this differential equation to be well 

;/ is of the same order as m, and is very small. Therefore, 



defined, / must be a function of 0(1). Thus, / = — , 

the potential (oc mf oc 1/A 2 ) is very small. 

Finally, note that Q\/{a 2 n 2 ) is a function of k 2 /(a 2 K 2 ) and H 2 /k 2 . Because of Eq. (|1U5|I for FIR modes, we guess 
that 0/j may be almost independent of k for r)k < Vk- This implies that the initial mode solution created at 77° of the 
FIR modes are almost independent of k. 

Since m is small, small variation 8rj corresponds to a large variation in the mode dependent conformal time Srjk ~ A<5?7 
in this region. Since the multiplication of potential and the time-interval, 



4>k 



A' 



(110) 



is small, the mode solution does not change much until r]^ after the mode generation at r\i. Then, the state at 77^ = 77^ 
is not much different from its ground state at rjj.. 




FIG. 3: Schematic plot of fi fc for IR modes. . 

Collecting all of the above arguments, we propose that the initial mode solutions, which are almost independent of 
k, vary only by a small amount during 77J? < 77 < 77^ ~ rj m . The evolution for 77^ > 77^ will be described by the same 
method as that of the UV modes with the matching condition at 77^ . The mode solution is given by Eq. I|75|l and its 
coefficients are given by Eq. |76{l. For these modes we have —krj^(^ —kr\ m ) < n(= —k m r\ m ). In the case of a FIR 
modes, we have \krj^ \ <C 1. Then, we can use the form l|77|l of the Bessel function to determine 



?r(t/+l/2) / k 



-1/2 



(111) 



where we have assumed T^f^T^) <C 1. 

Thus the power spectrum for the FIR modes is given by 



P F i R {k) 



n{v + l/2) 2 



32tt 2 is 2 a 2 (?7 m )M2fi fe (77 m ) 



(112) 



Since 77^ ~ r) m , we have strong (proportional to k 3 ) blue spectrum for this FIR modes. Since [i — v is a small positive 
number, the size of the spectrum slowly decreases in 77^. It was argued that, in relation to the pre-big-bang scenario, 
this kinds of cutoff of power spectrum in the low frequency region can explain the low CMB quadrupole moments |13| . 

For the spectrum to be continuous, the red spectrum in the UV region must be continuously deformed to the blue 
spectrum in the FIR region. Therefore, the spectrum for a$K < k < k m must be continuously deformed from the 
weak red spectrum to flat and then to weak blue spectrum. 
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VI. SUMMARY AND DISCUSSION 

We have studied the efffects of the k— deformation of Robertson- Walker space on the evolution of metric fluctu- 
ations in expanding cosmological background. For a given noncommutative n— deformed inflationary universe the 
cosmological background is still described by the Einstein equation since the background fields only depend on one 
variable t so that the homogeneity and isotropy of the Robertson- Walker space are kept. The equation for linear 
fluctuations, however, are modified. We have shown that the modification takes the form of nonlocal interaction of 
the fluctuating field with itself and with the background. We have analyzed the system by perturbatively expanding 
the action up to the first order H 2 /k 2 . 

An important consequence of the space-time non-commutativity is that for each wave number k, there exists an 
earliest time rfl at which the fluctuating mode is created. The origin of the mode generation phenomena is a direct 
consequence of the k— deformation which introduces an upper bound of the comoving wave-number by Eqs. (|12|l and 
(|67|l . We assume that the fluctuation starts out with its vacuum amplitude at rfl since the number of excitations 
should be conserved during the creation process. Moreover, this condition restricts the physical frequency to be smaller 
than k. This condition determines the initial condition of given modes for a given initial time. The deformation also 
generates the correction terms proportional to H 2 /k 2 and k 2 /(K 2 a 2 ) to the frequency squared, which determines the 
time evolution of mode k. 

There are two main results for the corrections to the power spectrum of the metric fluctuation due to the deformation. 
The first is that the deformation alters both the time dependence and the momentum dependence of the power 
spectrum. Especially, in the case of a power law inflation, we have shown that the power spectrum slowly decreases 
in time as {rjk/rjo) 2 ^~ u \ where ji — v~ cti 2 is a very small but positive. This time dependence of the power spectrum 
is a new feature of the present approach in contrast to the result of the commutative space case and to the result of 
Ref. 01 ■ The fact that the power spectrum decreases in time in addition to the existence of the red shifted spectrum, 
can be used to resolve the gravitational hierarchy problem for the tensor mode fluctuation. Another consequence of 
the deformation is that the momentum dependence of the power spectrum for UV modes is also dependent on time 
as k 3 ~ 2 ", where v is given by Eq. (92). Note that 3 — 2v is a small positive number for rjf < rjfi at and is a small 
negative number for rjf > rjfiat- Therefore, there is a period (77 < rjfiat) of blue spectrum in the earlier time of the 
inflation and the spectrum becomes red after the time rjfiat-, which is determined by n. Maximal red shift occurs to 
the power spectrum as rjf — > 00 in which case ji = v. Another interesting effect appears in the power spectrum of 
infra-red modes. Although we cannot obtain the dynamics of infra-red modes explicitly, we have suggested a form of 
the power spectrum from the consistency requirements. The power spectrum of the far infra-red modes (k < aore), 
with ao the scale factor at t = 1/k, have a cutoff proportional to k 3 even though the explicit procedure needs much 
refining since the dynamics at early times are not known. One can use the existence of this cutoff to explain the low 
CMB quadrupole moment. For the spectrum to be continuous, this A; 3 type power may change as k increases. We 
know that for ultra-violet mode {k > k m , Eq. (89)) the spectrum becomes slightly red shifted for rj > rjfiat- Therefore, 
the spectrum may change from weak blue to flat spectrum for a^n < k < k m . 

Brandenberger and Ho |T(ij computed the effect of the stringy space-time uncertainty relation to the power spectrum 
of metric fluctuations for power law inflation in the Robertson- Walker space. It is interesting to compare their results 
with that of the present paper since we have started from a different commutation relation (|23[) . We start from 
the basic commutation relation l|23ll and construct the theory from the first principle, although we have to use the 
perturbation to compute the physical effects of the deformation. Let us consider the power spectra for time rjf > rjfiat- 
The spectra of Ref. [T(| changes from weak red for ultra-violet modes to weak blue for infra-red modes. In our case, the 
far infra-red modes behaves as k 3 . This is due to the fact that the far infra-red modes becomes almost independent of 
k for the very early times rj < rj m . Since the generation time for the ultra-violet modes is of similar form as in Ref. |10| , 
the behavior of the spectrum in our case is similar to theirs for UV modes. The spectrum for k > a^n changes from 
weak blue to weak red as k increases. A totally new phenomena due to the k— deformation, which is absent in the 
cases of Ref. 0] and the metric fluctuations in the commutative Robertson- Walker model, is the time-dependence of 
the spectra. We have shown that both the spectra of the ultra-violet and infra-red modes decrease slowly in time. 

We have used the perturbative approximation to obtain the generic feature of the physical effect of the 
k — deformation on the cosmological evolution. It would be interesting if one could develop some nonperturbative 
approximation methods to extract better information from the nonlocal theory described by the action l|29|l . 
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